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Abstract 

' A Lagrangian derivation of the Equations of Motion (EOM) for static spherically symmet- 

, ric metrics in F{R) modified gravity is presented. For a large class of metrics, our approach 

permits to reduce the EOM to a single equation and we show how it is possible to construct 
exact solutions in _F'(_R)-gravity. All known exact solutions are recovered. We also exhibit a 
new non trivial solution with non constant Ricci scalar. 



o . 

O^. 1 Introduction 

bJO, Recent observational data imply that the current expansion of the universe is accelerating. This is 

the so called Dark Energy issue. There exist several descriptions of this acceleration. Among them, 
^■f-^ , the simplest one is the introduction of small positive Cosmological Constant in the framework of 

^ ' General Relativity (GR), the so called A-CDM model. A generalization of this simple modification 

, of GR consists in considering modified gravitational theories, in which the action is described by 

■ a function F{R) of the Ricci scalar R (for example, see [T]-^]). Typically these modified models 

I admit the de Sitter (dS) space as a solution and the stability of this solution has been investigated 

• ■ in several places (see, for example [HllllElE]). Furthermore, viable F{R) models, namely models 

' able to pass the local gravitational GR tests as well as able to describe the inflation with dark 

energy in unified way, have been recently discussed [71 |8l |9l [10] . 

To begin with, we recall the Equations of Motion in vacuum for a general F(R) model 



X 



Gfiv = Rfiu — 2^9tJ.i' — . (1) 



■ Here, R^i, is the Ricci tensor and the 'modified gravity' tensor is given by 

C = [Id^Am) - RF'm + (v^v. - <?^.n)F'(i?)| . (2) 

The prime denotes derivative with respect to the curvature R, is the covariant derivative 
operator associated with g^^ and Dip = g^'^V ^jV u4) is the d'Alembertian. The trace of Eq.((T]) 
reads 

3DF'{R) + RF'{R)-2F{R) = 0, (3) 

which shows that there exists an additive scalar dynamical degree of freedom represented by F'{R), 
absent in GR. 

The issue to find exact static spherically symmetric (SSS) solutions different from the Schwarzschild- 
dS one of these modified gravitational models appears a formidable task, since also for a reasonable 
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F{R), the equations of motions are much more compUcated with respect of the ones in "vacuum", 
namely absence of matter, of General Relativity. For the specific choice a class of exact SSS 

solution has been presented in 

A general discussion on SSS solutions based on above equation of motions has been presented 
in [121 [T^ [HI I15| . where one can find further references. 

In this paper, we would like to present an alternative method for constructing SSS of a generic 
modified gravity models, and, as a result, we also exhibit new exact SSS solutions. 

The paper is organized as follows. In Sec. II, we present our Lagrangian derivation of the 
Equations of Motion for SSS metrics in F(i?)-modified gravity. In Sec. Ill, a large class solutions 
are investigated. In Sec. IV, we study the highly non trivial case of solutions with non constant 
scalar Ricci curvature. Finally, conclusions are given in Sec.V. 

We use units offc^ = c = h = 1 and denote the gravitational constant = SttGn = STr/Mp^ 
with the Planck mass of Mpl = Gj^^^^ = 1.2 x lO^^GeV. 



2 Lagrangian approach for static spherically symmetric vac- 
uum solutions 



To begin with, we recall the action for modified F(i?)-theories: 



(4) 



where g is the determinant of metric tensor g^^ and F(R) is a generic function of the Ricci scalar 
R. 

We shall look for static spherically symmetric (SSS) solutions of the type. 



ds^ = -B(r)c2°('-)di2 + 



'dn. 



(5) 



where dVl = r^{d9^ + sin^ 9d(p) and a(r) and B{r) are unknown functions of r. With this ansatz, 
the scalar curvature reads 



R = 



dr 



B{r) 



dr 



a{r)-2B{r) 



a (r) 



dr 
B{r) 



^Bir)^2B{r)^air) 



(6) 



By plugging this expression into the action ^ , one obtains a higher derivative Lagrangian theory. 
In order to work with a first derivatives Lagrangian system, we may use the method of Lagrangian 
multipliers, used for the FRW space-time in Refs.flOllTTllS]. In the static case we are dealing with, 
the method permits to consider as independent Lagrangian coordinates the scalar curvature R, and 
the quantities a and B, appearing in the spherically static symmetric ansatz. As a consequence, 
we will obtain two equations of motion where one of the unknown quantities a appears in a very 
simple way. The main difference with respect to the other general appraoches is that we do not 
directly make use of EOM 

By introducing the Lagrangian multipliers A and making use of Eq. ([5]), the action Q may be 
written 



1 

2^ 



dt / dr 



(e^f'^V^) |F(i?)-A|i?+|; 



dr 



B{r) 



dr 



a (r) 



-2B{r) 



dr 



B{r)f,a{r) 



— B(r) + 2B(r)— a(r)+4* 



B{r) 



rrB{r) 



(7) 
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Making the variation with respect to R, one gets 



\ = F'{R). 



(8) 



The prime denotes the derivative with respect to the curvature R. Thus, by substituting this value 
and by making an integration by part, the Lagrangian takes the form 



L{a, da/dr, B, dB/dr, R, dR/dr) 



" <^ - F'{R)R) + 2F'{R) 1 



dB{r) 
dr 



Bir 



dR^fdBir) 
dr 



+ 2B{r) 



da{r) 



dr dr 

Making the variation with respect to a, one gets the first equation of motion 
RF'{R) - F{R) ^ {l- B{r) ~r{dB{r)/dr)) 



F'iR) 
2B{r)F"{R) 
F'{R) 



(9) 



(10) 



d'^R /2 dB{r)/dr\ dR F"'{R) ( dR 
TP" ^ Vr ^ 2B{r) J d7 ^ F"{R) V d7 



= 0. 



The variation with respect to B{r) leads the second equation of motion 



da{r) /2 F"{R) dR\ F"{R)d^R F"'{R) fdR 
dr \r^ F'(R) dr J F'{R) dr^ \ dr 



= 0. 



(11) 



while by making the variation with respect to R, we recover the Eq.®. 

Once F{R) is given, together with equation ©, the above equations form a system of three 
differential equations in the three unknown quantities a{r), B{r) and R{r). We would like to note 
that one advantage of this approach is that the a does not appear in Ea. ljlOp and vice versa. In 
the next Sections, we will find exact solutions of the above system of differential equations related 
to specific choices for F{R). 

As first check of the formalism, let us look for the well known case of constant curvature 
solutions R ~ Ro- From Eq. (fTT|) we have a — Const and it is easy to show that the only solution 
of Eqs.® and ([TU)) is the Schwarzshild-De Sitter solution 



where C is a generic constant and 



with Rn = 4A. 



Bir)^l-^r' + ^ 
6 r 



F{Ro) 



A = 2 ( i?o 



F'iRo) 



(12) 



(13) 



3 Solution with constant a 



Now, let us consider the case of non constant Ricci curvature, but still with a — Const. From 
Eq.dni): 



Thus, 



F'{R) ^ar + b. 



(14) 
(15) 
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Here, a and b are two integration constants. If we give the explicit form of R, we may find r as a 
function of Ricci scalar and reconstruct F'{R) realizing such solution. The equation ^ leads to: 

Since {F"{R))dR/dr = dF\R)/dr = a and dF{R)/dr = F'{R)dR/dr, by multiplying Eq.(fTU| 
with respect to F'{R) and by deriving with respect to r, one has 

d'^B(r) [ h\ 2a , , , , 2h , , , , adBir) , , 

'- + 1 +ij(2^W-i) + i^(^W-i)-r-7r^ = 0- (17) 



dr^ 



r 



If a = 0, the solution of this equation is the Schwarzschild-de Sitter solution, which corresponds 
to F'{R) = b, namely GR plus a cosmological constant F{R) = bR + Fq. If 6 = 0, the general 
solution is 

Bir)^l(l-^ + C,r-). (18) 

Eq.llll) leads to 
so that 



nR) = a,lj^^. (20) 



We then can reconstruct F{R) and we get the model 



F{R)^-y/R + 6C2, (21) 

for which Eq. (fT5)) is a spherically symmetric static solution with a ~ 0. If we put C2 = 0, we 
recover a special case of Clifton-Barrow solution [TT]. If Ci — 0, we recover a special case reported 
in [HI. 

The most general solution of Eq. (IT7|) for 6 7^ is 

^, . f r^C2log{r)a^ C 2 log {b + ar) Sr^a^ rC2a^ log (r) 

= 1 b^ + b^ + b^ b^ (22) 

log {b + ar) ra Caa ^'a , 

+ b^ T+ 26^+'' ^^"3&;; + ^ 

where Ci and C2 are generic constants. For a — 0, one again obtains the Schwarzschild-De Sitter 
solution. The above solution is one of the main result of the paper, because it is the starting point 
of the reconstruction method, and it is compatible with the result obtained in Ref. |15) . Implicitly 
F{R) is determined by Eqs. ([T5|) . and (22). If a F(i?)-model realizes the metric 

ds^ = -B{r)df + —— + r^dQ^ , (23) 



the coefficient B{r) assumes the generic form of Eq.(| 

On the other hand, we have to note that, since Ea. (jl7l) has been obtained trough a derivation, 
the left side of Ea. ([TU)) evaluated on the solution ([2^ could be proportional to 1/F'{R), so that 
some other constraints on free parameters could be necessary. As an example, let us consider the 
simple case where 

b^ 

C2 = . (24) 

a 

Since b is adimensional, we assume b — 1 and Ea. ip^^ becomes 
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and the Ricci scalar reads 



By using Ea. ([T5|) one has 



R=\- 12Ci - 18a^ . (26) 



F{R) = R + 2ay/R + 18a^ + 12Ci . (27) 

Now, we find that the left side of Eq. (|10p is 3(3a^ + 2Ci)/F'(i?), so that we have to require 
Ci = — 3a^/2. As a result, the reconstruction gives 

F{R) ^ R + 2ay/R, (28) 

and this model admits a spherically symmetric solution with constant a parameter and 

+ , (29) 



2 V 3ar_ 
and Ricci curvature R = l/r^. 

4 Solution with non constant curvature 

In this Section, we will look for solutions with non constant a and non constant Ricci scalar. Thus, 
we make the ansatz 

« = ilog(^)\ (30) 

i? = (5r% (31) 

where ro,q,S and s are constant parameters. From Eg. pip , it follows that the modified gravity 
model which realizes this kind of solution is 

F{R) ^ k{R^ + ^i) , (32) 

where k and /x are constants and 7 is given by 



(2 + g + 4s)s± Vg2 + 20g + 4 
T = -^2 ■ (33) 

If we put s — —2, one has q — 2(7— 1)(27— l)/(2 — 7). In this case, Ea. ([TU)) has a simple solution 
by choosing /j, = and 5 = 67(7 — l)/(27^ — 27 — 1). Thus, one recovers the Clifton-Barrow 
solution |TT] 

^^""^ " (107- 472 - 7)(272 - 27- 1) " r(47^-i07+7)/(2-7) ) ' ^^4) 
which is consistent with Eq.®, compatibly with assumption (|31l) . 

However, there exists also another solution with q = 1 and s = — 1. It means, 7 = — 1 or 
7 = 3/2. Eq.® is solved by 

Here, Ci and C2 are arbitrary constants. On the other hand, it is easy to see that Eg. pTI)) is 
inconsistent if 7 = 3/2, but admits a solution for C2 = if 7 = — 1 and 



(36) 



, 2 

with /i negative. As a consequence, we may put /i = — |-, /i > and we have found that the 
model 



R 6 



admits the following SSS metric 



with 



The Ricci scalar becomes 



ds' = -i -) B{r)dt' + 

Jo J B{r) 



A f 7 c 



h r 



r^dn , 



(37) 

(38) 

(39) 
(40) 



5 Concluding remarks 

In this paper, a new method for investigating SSS solutions for a generic F{R) has been presented. 
It is our opinion that our method is quite general and offers some advantages with respect the 
other general methods presented in [T^l [131 Ull US • In fact in the important case of a constant (in 
this case, we may put a = 0), we obtain a general form of the metric via the explicit expression 
for the quantity B{r). 

The static solutions describe a black hole as soon as there exists a real positive solution of 
B{r}j) = 0. If this happens, there exists an event horizon and the related Hawking temperature 
reads 

This is well known result, and it can be justified in several ways, like the elimination of conical 
singularity in the associated Euclidean metric, or making use of the tunneling method introduced 
by Parikh and Wilczek [H [TH]. 

For example, for the model F{R) = ^/ R + 6C2 of Eq. (|2ip . since the positive solutions of 
B{rH) = read 



-1± VI + 4C1C2 
2C2 



(42) 



one has to require Ci > and C2 > (for the plus sign) or C2 < and Ci < — 1/(4C2) (for the 
sign minus). 

For the other model in Eg. ipS]) . we have to chose a — —g < 0, g > 0, thus F{R) ~ R — 2g^/R 
and we have 

^ff = |- ■ (43) 

With regard to the other class of solutions we have discussed, the ones with non trivial a{r), 
the situation is different. A discussion concerning the Clifton-Barrow class has been presented 
elsewhere [20J. Thus we consider only the new solution (j39|) corresponding to the model F{R) — 

k ^-^ — . In this case, if C > 0, there exists a positive root of 



1 - — — 



0. 



(44) 



since dB{r)/dr < 0, namely is non vanishing. Thus, as in de Sitter space-time, we may take the 
absolute value of the surface gravity for the expression of Hawking temperature. Furthermore, 
besides the usual Killing vector associated with the static space-time, on general ground, it is 
also possible to introduce the Kodama vector, with an associated surface gravity, the Hayward's 
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surface gravity. A detailed discussion about this issue can be found in [211 As far as the 

other thermodynamical quantities associated wit these black holes solutions, the entropy can be 
calculated by the Wald method ^24, ,25, i26j , and one has [3] 

Sh = ^F'{Rh). (45) 

In the last example, we may choose the dimensional constant k negative in order to deal with a 
positive entropy. 

Much more difficult is the issue associated with the energy or mass of these black hole solutions. 
With regard to this, we refer to the papers [?71 [5D] . 

Finally we note that the functional freedom in the choice of F{R) gives the possibility of its 
cosmological reconstruction as it has been reviewed in [T]. In fact, following the construction 
developed in this paper and BH reconstruction outlined in Ref. [31], we can generate more black 
holes in F(R) gravity via its reconstruction. 
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